Contractor renormalization group theory of the SU(A^) chains and ladders 
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Contractor renormalization group (CORE) method is applied to the SU(7V) chain and ladders in 
this paper. In our designed schemes, we show that these two classes of systems can return to their 
original form of Hamiltonian after CORE transformation. Successive iteration of the transformation 
leads to a fixed point so that the ground state energy and the energy gap to the ground state can 
be deduced. The result of SU(A'^) chain is compared with the one by Bethe ansatz method. The 
transformation on spin- 1/2 ladders gives a finite gap in the excited energy spectra to the ground 
state in an intuitive way. The application to SU(3) ladders is also discussed. 
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The contractor renormalization (CORE) group 
method combines the contraction and cluster expansion 
techniques with the real space renormalization group ap- 
proach to solve the electron and spin lattice problemsi It 
was first applied to spin- 1/2 Hcisenberg chain and (1-1-1)- 
dimcnsional Ising model, and later to the frustrated 
antiferromagnets and the Haldane conjecture. The 
results are satisfactory and encouragingiiiS Since then 
the method has been applied to investigate low energy 
physics in many strongly correlated systems i^i^i^iSiL^ 
In this paper, we are concerned with a class of models 
showing that CORE is at its critical point, which means 
that the symmetry of the system is restored or the 
same form of Hamiltonian is reconstructed after the 
CORE transformation, just like the spin-1/2 Heisenberg 
chaini Undoubtedly, this method is not limited to such 
a kind of systems. Though in many systems the original 
Hamiltonian cannot be recovered, the lower energy 
physics are retained and studied successively after the 
truncation and transformation. 

We consider the SU(Af) chain and ladders in this work. 
We show that the same form of the Hamiltonian is recov- 
ered after dividing adequately the lattice into blocks and 
defining a truncation scheme, so that the CORE algo- 
rithm can be done recursively. The range-2 result for 
SU(iV) chain is compared with the Bethe ansatz solu- 
tion by SutherlandiS The comparison suggests that the 
CORE method can give good result especially for large 

and relatively larger blocks. We also present the re- 
sults of the real-space renormalization group (RG) the- 
ory, which usually agrees qualitatively with the one by 
range-2 CORE calculation!^ In many cases the latter can 
be regarded as a refined method on the former. The 
spin-1/2 ladders have attracted a lot of attention since 
the discovery of a finite spin gap in the 2-leg laddersiiSiii 
Another CORE scheme based on plaquette-dividing of 
the ladder had been applied to this systemiiiS Here we 
shall use a different scheme which shows a S" = 1 magnon 
gap in an intuitive way. Results up to range-3 are pre- 
sented. 

It was shown by Morningstar and Weinstein that the 
CORE scheme of three-site block partition and two-state 
truncation on SU(2) chain recover the original form of 



Hamiltonian. Then the resulting effective Hamiltonian 
can be solved iteratively and a quite satisfactory result 
can be obtainedi The recovery of the form of the Hamil- 
tonian owes highly to the SU(2) symmetry and an ade- 
quate designed CORE scheme. As a generalization, we 
found that their CORE scheme on SU(2) chain is a spec- 
iman picked out from a general CORE scheme on the 
SU(iV) chain. Though the SU(iV) chain had been ex- 
actly solved by Bethe ansatz method long time agOf^ it 
is still instructive to see how CORE works in the system. 

Let us start with a one-dimensional SU(iV) chain in 
terms of the exchange operator, H = Pj.j+i- Here 
we limit our discussion to the antiferromagnetic case by 
setting J = 1 > 0. For a SU(A'') system each site j has 
N quantum states |j, a) with (a = 1, 2, • ■ • , N). The ex- 
change operator Pj,j+i swaps two states on sites j and 
j + 1, i.e. Pj^j+i\j,a;j + 1,(3) ^ \j,(3;j + l,a). Usually 
Pj,j+i can be expressed in terms of the SU(iV) genera- 
tors as Pjj+i = Y.apJpij)Ja{3 + 1), where the opera- 
tors J'^ii) satisfy the SU(A^) algebra Jp{j), Ji^iji) = 

^ii' {^uJ/^U) ^ ^pJuiJ)j- Alternatively, Pjj+i can also 
be expressed by spin operators^^ii^ Many spin systems 
as well as spin-obital systems concerning SU(A'^) symme- 
try have been studied cxtensivelyii^iiSiiLi2iiS 

In the CORE scheme, the first step is to divide the 
original chain into a chain of blocks and retain adequate 
number of energy levels in each block. We found two 
obvious schemes to be readily applied to this system: 
one is {N ~ l)-site block patition with A^-state trun- 
cation (scheme A), and the other is {N + l)-site block 
patition with A^-state truncation (scheme B). The treat- 
ment on SU(2) case in Refi obviously falls into the 
scheme B with Pj,j+i = 2Sj ■ Sj+i + 1/2 when = 2. 
We will see the scheme B gives better results than the 
scheme A. The existence of the two schemes can be un- 
derstood from the single column Young tableaux with 
(A^ - 1) or (A^ + 1) boxes. In fact the SU(A^) model 
on both (A^ — l)-sitc block and (A^ -|- l)-site block have 
one unique A^-dimcnsional ground state space. We de- 
note the truncated space for a single block by = 
{\(t>j,i) , l'/'i,2) , • • ■ , |0j,7v)}- Then in the range-2 CORE 
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TABLE I: The coefficients and in Eq. (Q 
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calculation, we should retain appropriate N'^ low levels 
from the exact diagonalization of two blocks. All the re- 
tained low levels should have nonzero projection to the 
product space $j (8) ^j+i: so the eligible levels are not 
always the lowest ones. Fortunately this job is easy to be 
done due to the SU(7V) symmetry. The range-2 CORE 
calculation leads to the effective Hamiltonian, 



NT 



(1) 



where the sign =p corresponds to the two schemes A (— ) 
and B (+), Pj,j+i is a "new" renormalized exchange op- 
erator connecting blocks j and j + 1 after each block 
"contracts" to a single site. The coefficients Cip and Kzf: 
are listed in Table It can be confirmed the range- 
3 Hamiltonian will include another operator Pjj+2 and 
the range-4 Hamiltonian will include more operators like 



P, 



1 P.i-J+iP', 



3.j+l-^j+2J+3 



].j+2 Pj+l,j+3 1 P].j+3Pj+l,j+2 ■ 



Here we only give the range-2 results since higher range 
calculation will not change the physics. For the scheme 
A and B, we give the results for N = 3,4,5 and 2,3,4 
respectively. 

Successive application of CORE on Eq. ^ will lead 
the running coupling approaching a gapless fixed point. 
And no phase transition is observed. The ground energy 
is read out as 
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(2) 



where the sign =p corresponds to the two schemes. Fig. 1 
shows that the result of the range-2 CORE of the scheme 
B agrees quite well with the one by Bethe ansatz method. 
The numerical error can be reduced by higher range cal- 
culation. The range-4 result for TV = 2 by Weistein shows 
the error is reduced to — 0.0025i^ 

In fact the traditional RG gives an effecive Hamiltonian 
having the same form of Eq. . It can produce results 
consistent with CORE though not so goodi^ The two 
schemes above are still applicable and the corresponding 
coefficients can be found in Table The advantages of 
CORE are obvious. In many cases one can design more 
flexible schemes in CORE while selecting basic blocks and 



FIG. 1; (Color online) The ground energy of SU(A'') chain. 
The scheme B of CORE gives better results and the numer- 
ical errors are about -0.0106, -0.0006, 0.0021 for iV = 2, 3, 4 
respectively compared to the results by Bethe ansata-. 



truncating at low levels i^^'^'SiS A more careful analy- 
sis shows that RG based on {N — l)-sitc block partition 
scheme (scheme A) can give an effective Hamiltonian for 
general N, 
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(3) 

which exhibits a ground energy coinciding with the one 
by Bethe ansatz method at large A^, = —N{N — 

2)/{N^-N + l)^-^ -1. 

The 2-leg spin- 1/2 ladders aroussed a lot of attention 
when a finite spin gap was observed^ A simple picture 
says that the ground state is a product state with the 
spins on each rung forming a spin singlet. Then the low- 
est energy excitation is a 1 magnon. Here we show 
that our scheme of CORE produces exactly the same 
picture and refined results can be achieved following the 
CORE algorithm. We start from the Hamiltonian, 
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where the indices A and B refer to the two rails of the 
ladders, a ~ Jmng/ Jmii is the ratio between the rung and 
rails couplings, and we have set Jraii = 1- 

Our first step is to divide the ladder into triads along 
the rail direction (Fig. 2(a)). The problem on the rail di- 
rection is just the SU(2) chain that had been solved. De- 
tailed calculation shows that the effective interaction be- 
tween the two blocks along the rung can also recover the 
Heisenberg interaction. Thus "new" ladders with renor- 
mahzed couplings can be obtained. The second step is 
to parse out the effective block-block interactions from 
all possible configurations of connected blocks. As de- 
fined by Morningstar and Weinstein, r connected blocks 
contain range-r' interactions with r' = 0, 1, • ■ • , r (r' = 




FIG. 2: (Color online) (a) 2-leg ladders. The basic blocks 
are triads along the rail direction. The fixed point is a chain 
of decoupled dimers. (b) The unsymmetric configuration of 
blocks involved in range-3 CORE calculation. 

corresponds to the constant term as in Eq. Q). To parse 
out all range interactions the exact diagonalization is em- 
ployed on the connected blocks. We present range-2 and 
range-3 results here. It is notable that the range-3 blocks 
should include a configuration in Fig. 2(b). This unsym- 
metric configuration may make the iteration procedure 
more troublesome. 

The range-2 CORE result simply regains the original 
form of Hamiltonian except for a constant term, 

H^'^ = (S/.S/+i+Sf Sf+i)+A(a) Sf-Sf], 

3 

(5) 

where 5 = 0.491582, C{a) and A(a) vary with a. The 
iteration on the range-2 effective Hamiltonian is always 
applicable because the retained four low levels are al- 
ways one spin singlet and three spin triplets just like 
the SU(2) chain case. After n steps of iteration on 
Eq. Q we will get running coupling terms as /i„ = 
<5„(S/ • + Sf • Sj'+i) + A„(a) §/ • Sf , where the 
coefficients are determined recursively, 5n = <5" , A„ {a) = 

<5"-i A(%#), ,A2(a) = 5 A(AM),Ai(a) = 

A(a),Ao(a) — a. So the rail coupling approaches zero 
(5" — > as n — > oo, while the rung coupling goes to 
a fixed value A„_^oo(Q!) ^ for a > (we observed 
that A„_,oo(q;) — > only when a = 0, which is in agree- 
ment with the conclusion drawn by DMRG(i2*2& and ex- 
act diagonalizationSi) . So the system flows to a fixed 
point exhibiting dimer covering on each rung of the lad- 
der. The spin gap is read out as As(a) = A„_,oo(a). The 
ground energy Eq is obtained by cumulating the constant 
term. Table Hll gives an example of iterations procedure 
for a = 1. 

The range-3 CORE result at the first run of iteration 
contains the next-nearest-neighbour interactions, 

+A(a) ^ ■ Sf + r!(a) (S/ • Sf+i + Sf • S/+i) 

+7(S/.S/+2 + Sf ■Sf+2)], (6) 

where C(a), b{a) and A (a) are different from the ones in 
Eq. (jSJ, 7 = 0.033975. 7 will vary with a in the suc- 
cessive iterations , 7„ = 7„(a), • ■ • ,71 = 7. After n-step 
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TABLE II: An example of the range-2 CORE iteration pro- 
cedure at a = 1. 




FIG. 3: (Color online) The ground energy and the gap for the 
spin- 1/2 2-leg ladder. The log- log plot shows that CORE and 
RG give correct gap in strong coupling limit q ^ 00. Data 
by other methods are adapted from Ref^i*— *— . 



iterations, we found that the only nonvanishing coupling 
is still the interaction along the rung A„^oo(q!) 7^ 0, so 
the physical picture obtained by the range-2 CORE does 
not change, i.e. the ground energy and the spin gap are 
produced in the same way. 

As we noted above, the unsymmetric configuration of 
blocks in Fig. 2(b) brings some troubles to the range-3 
CORE iteration. Unlike the SU(2) chain, the desired low 
levels may not always stay at the lowest positions during 
the iterations. And sometimes it is hard to select out the 
eligible set of levels from several possible candidates since 
each of them will lead to a recovered SU(2) symmetry. 
So different iteration procedures with different results are 
inevitable. When these situations take place, we resort to 
the principle: retaining the iteration procedure that gives 
the lowest energy^ although in our observations the val- 
ues of the results only have small difference. The range-2 
and range-3 CORE results for the ground state energy 
and the spin gap are illustrated in Fig. 3. For a com- 
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FIG. 4: (Color online) The ground energy of the 2-leg SU(3) 
ladders. The insets (a),(b) and (c) show the schemes used in 
the calculation . 



parison, data by other mcthodsiiiSiS^iS^ arc presented 
together. The ground energy agrees well with those by 
other methods in the whole range of interchain coupling 
a. This means that CORE algrithm can successively cap- 
ture the low energy physics of the system. The gap has 
relatively larger deviation at intermediate values of a. 
Nevertheless the discrepancy can be remedied through 
higher range CORE calculation. The range-3 gap is a 
little zigzag. This may be due to the unsymmetric con- 
figuration of range-3 blocks in Fig. 2(b). It is noteworthy 
that RG gives a gap simply as Af~^ = a, which captures 
the correct behaviour of the gap at strong coupling limit 
a — > oopii 

We also apphed CORE to the 2-leg SU(3) ladders, H = 
Y^j [iPj,A:j+i,A + P].,B-j+i,B) + a Pj,A-j,B\- The applica- 
ble schemes are presented in Fig. 4(a), (b) and (c). No- 
tice that all blocks are equivalent and a 3-state truncation 
is made in each scheme. The scheme (a) should be valid 
when the rung interaction a is large enough. While for 



small a, the schemes (a) (a < 1.0) and (b) (a < 1.58) 
are appropriate and scheme (b) is better than (a). All 
the three schemes lead to the fixed point with zero gap. 
We see that scheme (a) will be mapped to a SU(3) chain, 
which had been solved previously and gives a zero gap. 
And after the first mapping we applied 4-site block par- 
tition scheme on the chain in the successive iteration 
steps to produce the ground energy in Fig. 4. While 
scheme (b) and (c) will return to a new 2-lcg SU(3) lad- 
ders, hn = 5n{Pj,A-j+l,A + Pj,B-J + l,B) + A„(a) Pj^A;j,B, 

but we observed that the running couplings of the rail 
direction (5„ and the rung direction A„ will go to in- 
finitesimals of the same order— as we push the iteration 
steps to infinity, n — > oo, so a gapless phase is also ob- 
tained, which agrees with the result of scheme (a). The 
SU(3) model on 4-leg ladders can be analyzed in similar 
schemes and a gapless result is also expected. The result 
is reminiscent of the SU(2) model on a chain and on a 
3-leg ladders, which are also gapless. But unfortunately 
the above schemes or their analogues are not applicable 
for the 2-leg SU(4) ladders, which exhibites plaquette 
singlet-multiplet excitationiiSiSiSi One may have to re- 
sort to other kind of schemes. 

In conclusion, we have studied the SU(A'^) chain and 
ladders by the CORE schemes. We have shown that the 
effective Hamiltonian in the appropriate CORE schemes 
can regain its original form such that it approaches a fixed 
point by iteration of the CORE schemes. The ground 
state energy and the lowest excitations can be deduced 
from the fixed point. The results show that the SU(-/V) 
chain and the 2-leg SU(3) ladders are gapless, while the 
2-lcg spin- 1/2 ladder exhibites gapped phase originated 
from the rung dimmerization. 
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